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ABSTRACT 

.  \ 

The  use  of  sensors  configured  Into  various  geometrical  shapes  called 
arrays  are  Important  In  extracting  signal  Information  In  many  applications. 
This  paper  explores  the  effectiveness  of  a  two-dimensional  or  planar  array  In 
extracting  signal  Information  from  digitized  data  when  additive  Gaussian  noise 
and  additive  Impulsive  Interference  are  present.  Based  on  the  likelihood 
ratio  approach  the  data  output  of  each  sensor  Is  first  passed  through  a 
zero-memory  device.  The  zero-memory  device  Is  either  a  nonlinearity  or 
linearity  depending  on  whether  the  Interference  Is  present  or  not  present, 
respectively.  Further  processing  utilizes  a  receiver  which  converts  the  data 
at  the  output  of  each  zero-memory  device  Into  frequency  components  using  a 
dlscrete-Fourler  transform.  ,  The  frequency  components  are  then  weighted  and 
sunned  over  all  spatial  sensors.  The  weights  are  chosen  by  maximizing  a 
defined  performance  measure.  It  Is  shown  that  this  optimization  procedure 
does  not  take  Into  account  the  debilitating  effects  of  Impulsive 
Interference.  In  order  to  treat  theoretically  the  Inclusion  of  a 
nonlinearity,  to  combat  Impulsive  Interference,  the  concept  of  an  Ideal 
nonlinearity  Is  Introduced. 
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I .  INTRODUCTION 


In  many  applications  an  array  of  sensors  Is  employed  to  extract  useful 
Information  from  a  signal  propagating  In  a  medium.  The  Information  to  be 
extracted  may  be  the  temporal,  spatial,  or  spectral  properties  of  the  signal. 
Recently,  the  probability  distribution  Information  of  received  data  were  also 
shown  to  be  of  Interest  [1].  It  seems  natural  to  assume  that  the  optimum 
processor  would  be  configured  based  on  the  desired  Information  to  be 
extracted.  One  approach  to  design  an  optimum  processor  would  be  to  specify,  a 
priori  ,  the  Information  to  be  extracted  from  an  expected  signal  and 
configure  a  receiver  accordingly.  If  the  signal  Is  embedded  In  noise  this 
complexity  must  also  be  taken  Into  account.  In  this  paper  an  array  processor 
Is  defined  and  then  optimized  by  choosing  certain  unspecified  parameters  by 
maximizing  a  performance  measure. 

The  array  of  Interest  Is  configured  as  a  rectangle.  Although  the  Initial 
mathematical  details  do  not  require  an  array  configuration  to  be  specified, 
the  performance  results  do.  Therefore,  the  array  configuration  Is  specified 
at  the  outset.  As  shown  In  figure  (1),  there  are  M  sensors  in  each  horizontal 
row  and  N  sensors  In  each  vertical  column.  In  the  paper,  all  mathematical 
relationships  will  be  with  respect  to  the  top  left  corner  sensor  and  cartesian 
coordinates  will  be  utilized.  Examples  will  be  given  for  a  plane  wave 
arriving  In  the  plane  of  the  array. 
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In  environments  where  the  non-Gausslan  noise,  limiting  signal  detection, 
is  impulsive  theroetlcal  predictions  suggest  that  detection  performance  can  be 
significantly  improved  by  incorporating  nonlinear  zero-memory  devices  as 
preprocessors  to  the  standard  receiver  optimum  for  6auss1an  processes. 

However,  In  order  to  derive  optimum  tractable  results  from  a  likelihood  ratio 
formulation  somewhat  restrictive  assumptions  are  needed.  Nevertheless,  the 
usefulness  of  nonlinear  preprocessors  are  clearly  Indicated  In  these  analyses 
[2,  3,  4,  5,  6,  7].  A  processor  was  discussed  In  references  [8]  which 
utilized  a  nonlinear  preprocessor  based  on  the  estimation  of  quantiles  [2]. 
This  paper  discusses  the  performance  of  a  two-dimensional  array  processor 
which  Incorporates  either  a  linearity  or  an  Ideal  nonlinearity  [9]  at  each 
sensor  as  a  preprocessor. 

The  concept  of  an  Ideal  nonlinearity  Is  Introduced  In  the  time-domain  to 
facilitate  deriving  the  theoretical  performance  of  the  resultant  receiver  in 
the  frequency-domain.  A  similar  problem  was  discussed  by  Martin  and  Thompson 
[10].  They  Introduced  a  technique  to  Improve  spectrum  estimation  of  a  signal 
embedded  In  Impulsive  noise.  But  they  did  not  theoretically  evaluate  its 
performance.  The  advantage  of  Introducing  an  Ideal  nonlinearity  Is  that  this 
concept  allows  a  theoretical  performance  prediction  to  be  derived. 

From  Its  nature  Impulsive  Interference  tends  to  produce  a  broad 
spectrum.  For  this  type  of  process  classical  filtering  methods  are 
Ineffective  If  the  signal  Is  contained  In  the  same  band  as  the  Impulsive 
Interference.  However,  based  on  the  principle  of  the  likelihood  ratio,  under 
the  assumption  of  Independent  and  Identically  distributed  samples,  the  optimum 
processor  to  combat  Impulsive  noise  Is  one  which  utilizes  nonlinearities  as 
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preprocessors  [2,  3,  4,  5,  6,  7].  In  this  paper  several  examples  are 
presented  which  clearly  show  the  benefit  of  using  nonlinearities  as 
preprocessors  over  classical  filtering  when  impulsive  interference  limits 
signal  detection.  The  impulsive  interference  is  modeled  based  on  measurements 
from  a  certain  physical  phenomenon  which  is  shown  later. 

Before  presenting  the  processing  required  for  a  two-dimensional  array, 
frequency-wavenumber  spectrum  estimation  for  a  line  array  will  be  discussed. 
This  will  reveal  some  advantages  of  a  two-dimensional  array  over  a  line  array 
and  also  justify  the  array  processor  to  be  defined  In  the  next  section. 

Consider  the  coordinate  system  shown  in  figure  (2).  Distributed  in  the 
plane  are  several  sensors.  Their  positions  are  given  by  the  vectors  r^,  J  * 

1,  2,  ....  N.  In  addition  a  plane  wave  Is  propagating  in  the  plane  at  a 
velocity  and  direction  described  by  the  vector  TQ.  This  vector  Is  called 
the  slowness  vector  [11,  12,  13,  14].  It  points  in  the  direction  of 
propagation  and  has  a  magnitude  equal  to  the  reciprocal  of  the  velocity  of 
propagation,  1e,  |  «Q  [  -  l/*o,  where  Is  the  velocity  of  propagation  of 
the  plane  wave.  The  wavenumber  Is  related  to  through  the  formula,  Kq  = 
2wf«T},  where  f  is  frequency. 

The  data  received  at  the  J-th  sensor  is  given  by 
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The  Information  needed  In  this  pap°“  Is  contained  In  the 
cross-correlation  function 


E[x(tr  x(t2,  J2)]  -  Ets^  +  vrJ  )  s(t2  +  c^.r,  )] 


-t2**0,(r31  -rJ2)] 


where  It  Is  assumed  that  the  sensor  output  Is  space  stationary  [12,  16],  E 
represents  expectation,  and  R$[  ]  Is  the  autocorrelation  function  of  the 
signal. 


If  the  sensors  are  located  only  In  the  x-dlrectlon  then, 

E[x(tv  x(t2,  J2)]  =  R$[t1  -  t2  -  d  s1ne/*o] 

where  d  Is  the  separation  of  the  sensors  and  0  Is  the  arrival  angle  of  the 
plane  wave  with  respect  to  the  array. 

In  this  paper  It  Is  assumed  that  the  output  of  each  sensor  Is  available 
In  sampled  form  and  the  uniform  sampling  Interval  Is  one.  In  addition  all 
time  delays  will  be  Integral  multiples  of  the  sampling  Interval. 


61ven  a  line  array,  l.e.,  sensors  located  In  one-dimension  only,  the 
frequency-wavenumber  spectrum  estimate  for  the  discrete  case  can  be  defined  as 
a  two-dimensional  discrete  Fourier  transform  (DFT)  as  follows  [12,  17,  18,  19, 
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N0-l 

A(x,  K)  =  OAjH0)  ^  X(1,  J)  e"JiK  e“jJK  (1) 

where,  j  =  ^1,  x  -  2wr/N0  Is  the  normalized  frequency  [11],  r  =  0,  1,  2, 

....  Nq-1 ,  and  K  Is  the  normalized  spatial  frequency.  Here,  the  coordinate 
system  Is  defined  with  the  line  array  along  the  x-axIs  and  the  first  sensor  at 
the  origin. 

The  normalized  spatial  frequency  Is  represented  as 

K  =  2*  K'/N, 

where  K’  Is  defined  as  the  discrete  wavenumber  component  with  respect  to  the 
array,  and 


K'/N  *  (r/NQ)  (d/c)  sin 

*  d„  s1ne„, 
r  n 

where  dp  Is  the  sensor  separation  measured  In  wavelengths  and  c  Is  the  speed 
of  propagation  of  the  wavefront.  Therefore  the  arrival  angle  Is  given  by 


and  the  following  relationship  must  hold 
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From  the  property  of  the  OFT  both  positive  and  negative  components,  ie, 
frequencies  or  spatial  frequencies,  can  be  defined.  But  negative  frequencies 
are  redundant.  However,  both  positive  and  negative  spatial  frequencies, 
associated  with  positive  frequencies  are  meaningful .  Therefore,  K1  will  take 
on  both  positive  and  negative  values  according  to  equation  (2). 

An  example  of  the  frequency-wavenumber  spectrum  is  shown  in  figure  (3). 
This  type  of  presentation  of  the  frequency-wavenumber  spectrum  has  been 
considered  in  reference  [21].  A  two-dimensional  fast  Fourier  transform  (FFT) 
was  employed  to  produce  the  plot  in  the  figure.  Specifically  a  128  point  FFT 
and  64  point  FFT  were  used  to  generate  the  frequencies  and  wavenumbers, 
respectively. 

One  disadvantage  of  the  line  array  Is  that  Its  beamwldth  at  end  fire 
(e  «  ±  90°)  Is  broader  than  Its  beamwldth  at  broadside  (e  =  0°)  [22].  This 
can  be  seen  In  figure  (3)  by  counting  the  number  of  wavenumber  bins  for  each 
30°  Increment  starting  at  broadside.  There  are  16  bins  from  broadside  to  30° 
12  bins  from  30"  to  60°,  and  only  4  bins  from  60°  to  end  fire.  The  same 
relationship  holds  for  the  negative  wavenumbers.  In  addition,  as  frequency  i 
reduced  the  number  of  bins  between  the  two  end  fire  positions  also  reduces 
according  to  equation  (2). 

In  the  next  section  the  array  processor  for  a  two-dimensional  array  will 
be  defined.  It  will  be  a  generalization  of  the  discrete  frequency-wavenumber 
spectrum  estimate  discussed  above.  However,  the  defined  processor  is  also 
based  on  the  work  developed  In  the  reference  [23,  24,  25,  26]. 
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II.  ARRAY  PROCESSOR 


A  two-dimensional  array  of  omnidirectional  sensors  Is  shown  In  figure  1. 

The  elements  are  designated  as  XLJ,  L  ■=  1,  2 . N,  J  ■  1,  2,  ....  M  and 

the  spacing  between  elements  Is  denoted  by  d  In  the  horizontal  direction  and 
denoted  by  1  is  the  vertical  direction. 

As  shown  In  figure  1,  a  plane  wave  propagating  In  the  plane  of  the  array. 
Is  coincident  with  the  first  sensor  and  Its  wavefront  makes  an  angle  of  © 
degrees  with  respect  to  the  horizontal  direction.  The  time  It  takes  the 
wavefront  to  reach  sensor  XLJ  with  respect  to  X^  can  be  computed  from  the 
equation 

*LJ  ”  "  T)005®/0  +  d(J  _  1)sine/c  (2 

\ 

One  advantage  of  a  two-dimensional  array  over  the  line  array  is 
Immediate.  For  example,  let  L  *  1,  then 

ijj  -  d(J  -  1 )sine/c 

If  ©  *  0°  or  180°  the  same  value  is  obtained  namely,  I.  =  0.  In 

I  J 

general  there  Is  an  ambiguity  associated  with  predicting  the  direction  of 
arrival,  e,  of  a  plane  wave  with  a  line  array.  For  this  case, 

*  d(J  -  l)s1n©/c,  the  ambiguous  angles  are  denoted  by  the  set 
(©,  180°  -  e).  For  the  other  case,  «Ll,  the  ambiguous  angles  are  denoted  by 


the  set  (e,  -0) . 
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These  ambiguous  angles  are  completely  resolved  using  equation  (3)  because 
an  unique  set  of  time  delays  are  obtained  for  all  arrival  angles  within  the 
plane  of  the  array.  However,  an  ambiguity  will  again  arise  In  predicting  the 

direction  of  arrival  If  the  plane  wave  is  propagating  out  of  the  plane,  say  at 

some  angle  «.  Where,  when  t  =  90°  the  plane  wave  is  propagating  within  the 
plane  of  the  array.  Then  the  right  side  of  equation  (3)  will  be  multiplied  by 

sine.  As  is  well  known,  this  amlblgulty  can  also  be  resolved  by  utilizing  a 

3-dlmenslonal  array.  But  this  generalization  will  not  be  needed  in  this 
paper.  However,  this  ambiguity  may  also  be  resolved  by  utilizing  additional 
Information  concerning  the  possible  propagating  direction  [27]. 

As  stated  In  the  last  section  the  data  are  assumed  in  digital  form.  Let 
the  signal  and  noise  be  wide  sense  stationary,  mutually  independent  and 
zero-mean  processes.  The  Interference  will  be  a  zero-mean  deterministic 
process.  The  data  at  the  output  of  each  sensor  are  modeled  in  the  form 

x(  1 ,  L.  J)  =*  s(  1 ,  L,  J)  +  n(1 ,  L,  J)  +  1(1,  L,  J) 

1  =0,  1 . Nq-1 

L  =  1 ,  2,  ....  N 
J  =  1,  2 . M 

where  n(1,  L,  J)  Is  spatially  Independent  Gaussian  noise,  1(1,  L,  J)  is 
Impulsive  Interference  propagating  from  a  specific  direction  which  may 
coincide  with  the  arrival  direction  of  the  signal.  The  interference  is 
modeled  as  an  impulse  function  of  the  form 
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W— 1 

1(1,  L,  J)  =  X  Ak{«[i  -  (ig+  kr  +ylj)]  -  4[i  -  (ig+l  +  kr  +  TLJ)]}(4) 

where  w  Is  the  number  of  Impulse  functions  occurring  in  the  data  interval. 
Is  the  separation  of  Impulse  functions,  and  ylj  represents  the  appropriate 
time  delays  for  the  Interference.  The  symbol  6[  ]  is  the  Kronecker  delta 
function,  and  ig  is  the  temporal  location  of  the  first  Impulse. 

The  information  needed  is  contained  in  the  cross-correlation  function, 


E[x(1r  lr  x(12,  L2,  J2)] 


W  '2  -  <  V, 


♦  Rn<1,  -  12)  »tL, ^  -  L2  J2] 


w-1  w-1  (  , 

E  Z  \  \  «i, 

k,-o  k2=o  K2  r 


<vv*v,n 


+  1  ♦  k-|T  ♦  T| 


{‘t’2  -  (,g  +  k2r  +  -  <[,g  +  1  +  V  *  'rl,J,)]}|<5) 


2  2 


The  array  processor  is  defined  in  the  following  way  [23,  24,  25,  26], 


A(x) 


N 

E 


L=1 


H  ^  ^ 

E  A.  ,(X)  (1/JnJ  E  X(1,  L,  J)  e~jx1 
J=1  u  v  0  1=0 


(6) 
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+  The  symbol  means  that  the  cross  terms  are  zero. 
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where  the  symbols  are  as  defined  previously.  Here  we  utilize  complex  weights 
Alj(x)  which  differ  from  the  frequency-wavenumber  spectrum  estimate 
discussed  In  the  Introduction.  The  weights  will  be  chosen  to  maximize  a 
performance  measure.  But,  clearly  they  should  Include  the  proper  phase  shifts 
needed  to  steer  the  array  toward  a  signal  source.  If  the  signal  Is  not  a 
plane  wave  equation  (3)  will  not  give  the  proper  time  delays.  In  this  case  a 
search  for  the  appropriate  time  delays  can  be  mad.e  by  choosing  the  weights 
which  maximize  a  performance  measure.  This  procedure  Is  sometimes  called 
focusing  the  array  [22]. 

Later  we  will  Include  a  nonlinearity  In  equation  (5)  to  combate  Impulsive 
Interference. 

For  this  paper  only  the  expected  value  of  equation  (5)  will  be 
considered.  The  expectation  of  equation  (5)  reduces  to 


E[A(X)]  = 

h  h  4  J2  '  '  L2J2  ( 

t— '  -j^1 1  jXlo  ) 

O/N0)  JL  iL  E[X(1r  V  V  X<V  l2,  J2)]  e  e  >  (7 

1 1  1  o 


'1  2 


where  *  represents  complex  conjugate. 


The  term  In  the  bracket  represents  the  cross-power  spectral  density. 

This  will  be  derived  first  for  the  three  components  of  equation  (5).  They  are 
as  follows: 
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1}  Signal, 


O/N  ) 


'2  'v,  v* 


jxlj  jxl2 


RcM,  -  -  <  *.  i  -  ,  >1  •  '  • 


From  the  transformation  of  variables,  1?  ■  12,  v  =  ^  -  1?  and 

D  =  t,  ,  -  ,  we  obtain, 

L1J1  L2J2 


V1 


(1  -M>/N0)  Rs(*  -  0)  e^XD 


No-l-0 


«  e 


Jx0  S  (1'  -  |*  +  D|/N0)  R$(*)  e~jx* 


— Jx  1.  .  _Jx  li  i 

-1x0  /’v  L1J1  l2j2 

e  S(x) - *S(x)  e  1  1  e  1  i 

Hq  ->« 


(8) 


where  S(x)  Is  an  estimate  of  the  spectrum  of  the  signal  and  converges  to  the 
true  spectrum  as  N  -»*  [28,  29].  Therefore,  the  frequency  properties  have 

O  » 

been  separated  from  the  spatial  properties  In  the  limit  as  NQ  ■»«  . 
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(1/N0)  Z  Z  E[n(  1  ] ,  Lr  n(i2,  L2,  J2] 


-jxi..  jxi, 
e  e 


-  l2j2] 

0 


(9) 


where  N(x)  Is  the  spectrum  of  the  noise. 


3)  Interference, 


o'V  Z  Z 


-jxi,  jxi 

E[I(1r  Lr  ^ )  I(l2,  L2,  J2)]  e  e  2 


*  I(X) 


-jXy 

e 


LlJ  , 


jXy 


L2J2 


(10) 


where 


I(X)  -  (2  A2  w2/m0)  t,  -  cos(x)] 

and  we  let  A^  *  A.  I(x)  is  the  power  spectrum  of  impulsive  interference. 

It  depends  on  the  amplitude  A  and  the  number  of  Impulse  functions  occurring 
over  the  Interval.  Since  the  impulse  function  was  defined  to  have  zero  mean 
value,  the  term  [1  -  cos(X)]  is  present.  At  x  -  0,  I(x)  *  0  due  to  this  term. 
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The  other  term 


s1n(w\  t/2) 
wsin(\T/2) 

has  maxima  and  minima  occurring  over  the  frequency  band.  The  maxima  can  be 
found  from  the  equation 


(11) 


s1n(x-t/2)  =  0,  which  has  solutions  at 

x  -  2*  Iq/t  ,  IQ  -  1,  ±1.  ±2,  etc. 

When  the 'denominator  Is  zero  In  equation  (11)  so  Is  the  numerator.  The 
value  of  equation  (11)  at  the  solutions  can  be  found  using  L'  Hospital's 
rule.  At  the  solutions  equation  (11)  is  equal  to  one.  An  example  of  this 
function  will  be  given  later  In  the  paper.  Since  a  model  for  the  Impulsive 
Interference  Is  available,  it  may  be  possible  to  estimate  the  parameters  A,  W, 
and  rand  substitute  them  into  1(A)  to  obtain  an  estimate  of  the  spectrum  of 
the  Impulsive  Interference  directly.  But  It  will  be  clear  from  the  examples 
that  this  can't  be  done  without  error.  However,  this  paper  considers  the 
performance  of  using  nonllnearltles  only. 

Using  the  above  results  equation  (7)  reduces  to 
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(X)  A* 
l2 


-jX4 

•.(X)  e 
J2 


L1J1 


jxl 

e 


L2J2 


s  ,  (X)  A.  ,  (\) 

Lri  l2j2 


-JXy 

e 


L1J1 


jXY 

e 


L2J2 


The  weights  can  now  be  chosen  to  maximize  a  performance  measure.  But 
before  this  Is  done  let. 


jxl 

AU(X)  -  e 


LJ 


Then, 


E[  A (X) ]  =  S(X)  (NM)2  +  N(X)  NH 


”jX(  Y 
e 


L1J1 


)  JX(Y 
e 


L2J2 


l2j 


2 

In  the  above  equation  the  signal  Is  amplified  by  (NM)  and  the  noise  Is 
amplified  by  NM.  Whereas,  the  amplification  of  I(\)  depends  on  its  direction 
of  arrival.  If  the  Interference  Is  coming  from  the  same  direction  as  the 
signal,  then 


•/■ft 
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\  J1  'Vl  •  V2  'V2 


and 


E[ A (X) ]  =  [S(X)  +  I(X) ]  (NH)  +  N(X)  NH 


Therefore,  when  the  Interference  Is  coming  from  the  same  direction  as  the 
signal  the  array  passes  both  equally. 

These  Ideas  can  be  stated  more  precisely  if  we  define  a  performance 
measure. 

Let  the  performance  measure  (PH)  be  defined  as  follows: 


E[  A(X)|H  ]  -  E[  A (X)  J  H  ] 

PH - EiAwrv —  <13 

where, 

H-| :  x(  1 ,  L,  J)  -  s ( 1 ,  L,  J)  -  n(1 ,  L,  J)  +  1(1,  L,  J) 

Ho:  x(1,  L,  J)  -  n(1,  L,  J)  +  1(1,  L,  J) 

Rewriting  equation  (13)  In  a  form  convenient  for  optimization  [25,  26]  we 


obtain 
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where  the  term  In  the  bracket  Is  the  gain  of  the  array  and  Its  components  are 
defined  as  follows, 


(A-j<j(X),  A-^X) , 


Anm(X)) 


*t  jUn  ^S12 
V  ■  (e  ",  e 


u*T  -  (eJXT",  eJXT'2 . .JXT|«) 


where  T  represents  transpose  and  I  is  an  (NM  x  NM)  Identity  matrix. 


As  Is  well  known  [25,  26,  30]  the  optimum  gain  Is  given  by 


*T  -1 

G  »  V  Q  V 


and  the  optimum  weight  vector  Is 


TM  No.  851054 


_ U\) 

N(X)  +  I(X)  NH 


Therefore,  the  optimum  weight  vector  reduces  to 


A  =  u  _ _ im _ 

-  1  N(X)  +  I(X)  NH 


V 


and 


e  ,  uy _ HX) -  y*T  *T  w 

\  m  N(X)  +  I(X)  NH  -  -  -  - 

*T  *T 

Before  Interpreting  the  function  V  u  u  V,  let  It  equal  NH.  Then 
If  I(x)  NH  »  N(x) ,  G0  ~  NH  -  1.  On  the  other  hand  If  V*T  u  u*T  V  ~ 

(NH)2,  then  GQ  ~  0. 

Therefore,  the  optimum  gain  depends  on  the  value  of  the  function 
*T 

V  u  u  V.  This  function  can  be  Interpreted  as  the  array  beam  pattern  (BP) 
since 


N 

BP  *  V*T  u  u*T  V  -|(1/NH)  £ 


£ 


e-JX  ( B— 0)  2 


(14) 
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where, 

B  *  J(L-l)  coseB/c  +  d(J-l)  s1neB/c 
0  -  4( L-1 )  cose/c  +  d(J-l)  slne/c 

and, 

B-D  *  J(L-l)  (coseQ  -  cose)/c  +  d(J-l)  (s1neB  -  s1ne)/c 

In  this  formulation  the  array  Is  steered  In  the  direction  eg  but  the 
signal  arrives  from  the  direction  e. 

From  these  results  equation  14  simplifies  to 


~$1nfN\»a/2cl12  ' 
Ns1n[Ua/2c]  j 


s1nfMxdb/2c' 

Ms1n[Xdb/2c] 


where. 


a  *  coseg  -  cose 
b  -  cosefl  -  sine 

and  the  parameters  a,  b  are  bounded  by 


-2  <  a  <  2  and  -2  <  b  <  2. 
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The  spatial  location  of  a  signal  coincides  with  the  maximum  value  of  BP. 
Therefore,  for  a  specific  signal  direction  e,  the  function  should  have  only 
one  maximum  value. 

The  possible  maxima  of  BP  can  be  found  from  the  solutions  of  the  equations 

s1n(\la/2c)  =  0  and  s1n(\db/2c)  =  0 

which  reduce  to,  respectively, 

a  *  I  v  /\  and  b  «  I.  X  /d 

o  w  o  w 

where  IQ  *  0,  ±1,  ±2,  etc,  and  xw  Is  the  wavelength, 
suppose,  e  *  0  and  »B  *  * 

Then,  a  *  -2  and  b  »  0.  Therefore,  If  the  array  Is  not  designed  properly  BP 
can  have  a  maximum  at  e  *  0  and  eB  *  *.  This  phenomenon  Is  called  a  grating 
lobe  [22].  It  can  be  prevented  If  I  <  x^/2. 

Similarly  at  e  -  */2  and  ©B  -  3/2*,  then  a  -  0  and  b  *  -2.  So  If 
d  <  x^/2  a  grating  lob  Is  also  prevented  from  occurring. 

These  design  parameters  J  and  d  <  \^/2  also  prevent  the 

optimum  weight  approach  from  cancelling  the  signal  due  to  grating  lobes. 
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III.  IDEAL  NONLINEARITY 

From  the  previous  section  the  complex  weights  ALJ(x)  can  be  chosen  to 
maximize  the  gain  of  the  array  processor.  But  for  interference  arriving  from 
the  same  direction  as  the  signal  all  that  can  be  done  is  to  allow  both  signal 
and  interference  to  pass.  Filtering  is  also  not  effective  against  impulsive 
Interference.  This  can  be  seen  from  the  relationship 

E[|X0(X)|2]  =  |H(X)|2  E[|X(X)|2] 

2 

where,  E[|X(x)|  ]  represents  the  spectrum  of  the  data.  Since  the 
Interference  Is  Impulsive,  E[jx(X)p],  has  a  broad  spectrum  and  a  filter 
H( x)  would  operate  on  the  signal  as  well  as  the  noise.  However,  from  the 
likelihood  ratio  approach  for  1.1. d.  data  the  correct  procedure  is  to 
preprocess  the  data  through  a  nonlinearity  [2,  3,  4,  5,  6,  7].  We  will  employ 
a  nonlinearity  even  though  the  data  is  not  independent.  Related  work  in  this 
area  can  be  found  in  references  [31,  32,  33,  34,  35]. 

In  order  to  derive  a  theoretical  relationship  for  the  array  processor  we 
define  an  ideal  nonlinearity  as  our  preprocessor. 

Definition:  Ideal  Nonlinearity 


y(i.L.J)  - 

o  , 

at 

1  ■  'g  ♦  «+  \l 

y(i.L.J)  - 

o  . 

at 

'■V'*"'  TLJ 

y(i.L.J)  - 

X(1,L 

,J) 

,  otherwise 
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This  definition  can  be  equivalently  expressed  in  the  form 


w-1 

y(I.L.J)  -  z(i.L.J)  -  £  z(1,L,J)  (<[1  -  (1  +  kT+  y  )] 

k*0  '  9  LJ 

+  «C1  -  1g  +  1  +  k-C+  tlj)  J } 


where. 


z(1,L.J)  =  s(1,L,J)  +  n(i ,L, J) 


As  in  the  previous  sections,  the  information  needed  is  contained  in  the 
cross-correlation  function 


E [y(  1 1 ,L-| ,J-j )  y(l2.L2,J2)] 


E[z(11,L1,J1)  z(12,L2,J2)] 


w-1 


-E[z(11,L1,J1)  Z(12,L2,  J2)]  -  dg  +  k*+  YL  j  )] 


+  4[11  -  (1g  +  1  +  k*  +  yl  j  )] 


w-1  / 

£ 
k*0  * 


-E[z(1rL1,J1)  z(12,L2,J2)]  «[12  -  (1g  +  kt  +  yl  j  )] 


+  <ti2-(ig  +  i  +  kT:+  YL  j  ) ] 
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+  E [ z(  i -| ,  L-| ,  J-| )  Z(  1 2 » 1-2 » >^2 )  -^ 

£  E  l*M,  -  (1- +  k,T  +  yl  j  >1  +  «C1,  -('a  ♦  1  +  kt  *  irL 

k.*0  k2=0  1  9  11  9 

*  * 

x  { 4 C 1 2  “  (1g  +  k2*^  *  V2)]  +  4[i2  "  (1g  +  k2^"+  Yl2J 

For  simplicity  we  shall  let  ylj  =  £LJ,  l.e.^the  Interference  Is  coming 
from  the  same  direction  as  the  signal.  In  the  following  analysis. 

The  cross-power  spectral  density  estimate  follows  as 


-jx^  jxi? 


-jx  «L  j  jx  i,L  j 

{s(\)  -  S2<3  (X,w)  +  S4(x,w)}  e  1  1  e  22 

+  N(X)  -  N0  0(X,W)  +  N,(X,w)  \  4[L,J,  -  L,J„]  (16) 


d/N0)  EE  E[yd 


v  i 


J-,)  y(i2.  4.  J2)l  e 
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s2  3(\,w)  =  (i/N0) 


Gs(13.  k,  X)  «  R$ 

6.(1,.  k,  1.  X)  = 

S  3 

S4(X.  w)  -  (1/N0) 


N2  3(X,w)  and  N4(x,w)  are 
respectively,  but  with  R 


u-l  V'  -  5l2J2 


E  E 


k=0  i3  =  -5LJ 
3  L2J2 


[Ss(13.  k,  X)  +  Gs(13,  k,  1,  X)] 


V’-'Vi 

i  5  _r  [6s(i3.  k,  x)  +  Gs(13,  k,  l,  x)] 

3  4'L1J1 


i3  -  +  k1^  e 


jU13  -  (ig  +•  kt")] 


jx[i  -(1  +  1  +  kt  )] 

R  [1  -  (1  +  1  +  kt)l  e  J  9 

s  g 


~jx(k-i  -  ko)tr 

RsC(k1  -  k2)r  ]  e 


+  Rs[(k1  -  k2)t  +  1] 


-jX(k!  k2)te-jX 


+ 


k2)^ 


1]  e 


jx( k  -  k2) £  i 


of  the  same  form  as  S2  3(X,w)  and  $4(X,w), 
)  replaced  by  Ru(  ). 
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The  terms  S2  3(x,w)  and  N2  3(x,w)  are  functions  of  l.^  and 

L2J2.  But  as  Nq  •**!  their  effect  can  be  neglected  just  as  was  done 

for  the  cross-spectrum  of  the  signal. 

With  this  assumption  the  power  spectrum  can  be  separated  from  the  spatial 
properties  of  the  array.  Therefore,  the  expected  value  of  the  array  processor 
with  the  Ideal  nonllnearltles  Included  can  be  expressed  as  follows. 


E[  Anl(x)] 


ZESI  \,VX) V 


(X) 


h  J1  L2  J2 


(l/No)  5Z  Z  Etydvlv^)  y(12.L2,J2)]  e 

h  1 2 


-jxi,  jxi2 


[S( X)  -  S23(X,w )  +  S4(X,w)]  ZEZE  AMn  (x)  al7j 


~jx  Z,  j  jx?i 

LlJl  u 


h  Ji  L2  J2 


+  [N(X)  -  N2j3(X,W)  +  N4(X,W)] 


(17) 


Several  examples  will  now  be  given  for  the  power  spectrum  of  a  signal  In 
additive  Gaussian  noise  and  additive  Impulsive  Interference  without  and  with 
the  Ideal  nonlinearity  Included. 
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IV.  ILLUSTRATIVE  EXAMPLES 

Figure  4  shows  two  Impulse  functions  occurring  In  a  sample  of  data.  The 
first  burst  In  the  figure  Is  due  to  a  physical  phenomenon  but  the  second  burst 
appears  to  be  a  reflection  since  It  is  180  degrees  out  of  phase  with  the  first 
burst. 

Our  model  for  Impulse  functions,  equation  (4),  can  clearly  represent  the 
impulses  of  figure  4.  However,  the  model  Is  more  general  since  other 
functions  can  also  be  represented.  For  the  following  examples  a  special  case 
will  be  considered.  Figure  5  shows  the  time  history  for  our  first  example. 

The  data  consist  of  two  additive  sinusoids  In  addlvltve  non-white  Gaussian 
noise  and  two  additive  Impulse  functions.  One  sinusoid  has  a  level  of  20  dB 
below  the  other.  The  impulse  functions  are  separated  by  25.0  msec  and  have  a 
level  of  4  volts. 

The  corresponslng  spectra  are  shown  In  figure  6.  The  top  plot  represents 
the  spectrum  of  signal  and  noise  before  the  Impulses  are  added.  Since  the 
noise  Is  non-white  the  lower  level  sinusoid  Is  detectable  at  a  frequency  where 
the  noise  level  has  diminished.  The  next  plot  In  the  middle  of  the  figure  is 
the  spectrum  with  the  additive  Impulses  included.  These  results  are  predicted 
from  equation  (10).  After  the  Impulses  are  extracted  by  employing  the  Ideal 
nonlinearity  the  resultant  data  have  the  spectrum  shown  In  the  bottom  plot  of 
figure  6.  These  results  can  be  predicted  from  the  theoretical  results  of  the 
previous  section. 
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Even  though  there  Is  a  significant  Improvement  from  employing  the  Ideal 
nonlinearity  In  the  spectrum  there's  still  a  residual  degradation.  With 
additional  processing  It  may  be  possible  to  reduce  the  residual  degradation 
still  further.  However,  since  our  results  are  theoretical  they  can  be  used  to 
compare  other  methods. 

The  next  example  Is  shown  In  figures  7  and  8.  The  time  history  of  the 
data  Is  Identical  except  that  the  Impulse  functions  are  separated  by 
50.0  msec.  Since  the  separation  of  Impulse  functions  Is  doubled  there  are 
twice  as  many  maxima  and  minima  occurring  over  the  band  as  shown  In  the  middle 
plot  of  figure  8. 
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SUMMARY 


The  ability  of  a  two-dimensional  array  with  Ideal  nonlinear  elements  to 
extract  spatial  and  frequency  Information  from  a  signal  which  may  be 
Interferred  with  by  Impulsive  Interference  was  discussed.  After  some 
preliminary  remarks  on  discrete  frequency-wavenumber  estimation  the  array 
processor  for  a  two-dimensional  array  was  defined.  It  consisted  of  a  set  of 
frequency  dependent  complex  weights.  It  was  shown  that  the  weights  could  be 
chosen  by  maximizing  a  performance  measure.  This  procedure  would  be 
advantageous  when  the  wave  could  not  be  assumed  to  be  a  plane  wave  and  when 
Interference  was  present  but  propagating  In  a  direction  other  than  the  signal 
direction. 

For  the  Impulsive  Interference  considered  It  was  shown  that  filtering  was 
Ineffective.  However,  from  the  results  based  on  the  principle  of  the 
likelihood  ratio,  nonllnearltles  were  effective  against  Impulsive 
Interference.  In  order  to  discuss  the  effect  of  nonllnearltles  on  the 
spectrum  the  concept  of  the  Ideal  nonlinearity  was  Introduced.  This  allowed 
the  spectrum  to  be  theoretically  derived  when  ideal  nonllnearltles  were 
Included.  Two  examples  were  given  to  Illustrate  the  theoretical  results. 
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Figure  1.  Two-Dimensional  Array 
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